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Abstract 

We obtain an expression for the active gravitational mass (Tolman) 
of a source of the 7 metric, just after its departure from hydrostatic 
equilibrium, on a time scale of the order of (or smaller than) the hy- 
drostatic time scale. It is shown that for very compact sources, even 
arbitrarily small departures from sphericity, produce significant de- 
creasing (increasing) in the values of active gravitational mass of col- 
lapsing (expanding) spheres, with respect to its value in equilibrium, 
enhancing thereby the stability of the system. 
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1 Introduction 



One of the most remarkable aspect of general relativity, is the very special 
status that this theory confers to the spherical symmetry. This situation is 
particularly well illustrated by the Israel theorem. 

Indeed, as is well known pQ, the only static and asymptotically-flat vac- 
uum space-time possessing a regular horizon is the Schwarzschild solution. 
For all the others Weyl exterior solutions [2|, the physical components of the 
Riemann tensor exhibit singularities at r = 2m. Therefore, it is intuitively 
clear that for high gravitational fields, the evolution of sources of Weyl space- 
time should drastically differ from the evolution of spherical sources [3]. It 
is important to keep in mind that the sharp difference in the behaviour of 
both types of sources (for very high gravitational fields) will exist indepen- 
dently on the magnitude of multipole moments (higher than monopole) of 
the Weyl source. This is so because, as the source approaches the horizon, 
any finite perturbation of the Schwarzschild space-time becomes fundamen- 
tally different from any Weyl solution, even when the latter is characterized 
by parameters whose values are arbitrarily close to those corresponding to 
the spherical symmetry. This point has been stressed before jlj, but usually 
it has been overlooked. 

Notwithstanding, spherical symmetry is a common assumption in the 
study of compact self-gravitating objects (white dwarfs, neutron stars, black 
holes ) , furthermore in the specific case of non-rotating black holes, spherical 
symmetry should be "absolute" , according to Israel theorem. Therefore it is 
pertinent to ask, how do small deviations from this assumption, related to 
any kind of perturbation (e.g. fluctuations of the stellar matter, external 
perturbations, etc), affect the dynamics of the system?. 

In a recent work [3] it was shown that for a non-spherical source (even in 
the case of slight deviations from spherical symmetry), the speed of entering 
the collapse regime decreases substantially, as compared with the exactly 
spherically symmetric case. Also, the sensitivity of the trajectories of test 
particles in the 7 spacetime, to small changes of 7, for orbits close to 2m, 
has been brought out [6 . 

It is the purpose of this work to study further the behaviour of axysym- 
metric sources for very high gravitational fields. This will allow us to put in 
evidence the role played by the non sphericity (however small) of the source, 
on the outcome of evolution. To do so we shall obtain an expression for the 
active gravitational mass of the source, immediately after its departure from 
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equilibrium. Here "immediately" means on a time scale of the order of (or 
smaller than) hydrostatic time scale -see section 4 for more details. 

As initial configurations, we shall consider two interior metrics. These 
ones were found some years ago by Stewart et al. [7 , following a prescription 
given by Hernandez 8J allowing to obtain sources of Weyl space time, from 
known spherically symmetric interior solutions. 

The configurations to be considered are sources of the so-called gamma 
metric (7-metric) jH], JU|. This metric, which is also known as Zipoy-Vorhees 
metric [TT], belongs to the family of Weyl's solutions, and is continuously 
linked to the Schwarzschild space-time through one of its parameters. The 
motivation for this choice is twofold. On one hand the exterior 7-metric cor- 
responds to a solution of the Laplace equation (in cylindrical coordinates) 
with the same singularity structure as the Schwarzschild solution (a line seg- 
ment 9j). In this sense the 7-metric appears as the "natural" generalization 
of Schwarzschild space-time to the axisymmetric case. On the other hand, 
the two interiors considered have reasonable physical properties and general- 
ize important and useful sources of the Schwarzschild space-time, namely the 
interior Schwarzschild solution (homogeneous density) and the Adler solution 
|12j . All this having been said, we would like to emphasize that our main 
goal here is not to model the behaviour of a specific type of compact object, 
but to illustrate the effects of slight deviations from spherical symmetry, on 
the source. 

As we shall see, the value of active gravitational mass of the source eval- 
uated just after the departure from equilibrium, will be seriously affected by 
slight deviations from spherical symmetry. 

On the other hand, we are well aware of the fact that the 7 metric is 
not the only possible description for the exterior of a compact objetc and, 
of course, the two equations of state considered here do not exhaust the 
list of possible candidates for the equation of state of the stellar matter. 
However, in view of the properties of the 7 metric and the two equations 
of state considered here, mentioned above, it is fair to say that the case for 
the relevance of small deviations from spherical symmetry in the dynamics 
of compact objects has been further strengthened. This suggests that any 
conclusion on the structure and evolution of a compact object, derived on 
the assumption of spherical symmetry should be carefully checked against 
small deviations from that assumption, whenever the boundary of the source 
is close to the horizon. 
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2 The Space-time 



2.1 The exterior space-time 

As it has been mentioned above, our initial matter configuration is the source 
of an axially symmetric and static space-time (7- metric). In cylindrical co- 
ordinates, static axisymmetric solutions to Einstein equations are given by 
the Weyl metric [2] 



ds 2 



e 2X dt 2 - e 



with 
and 



e 2 ^ {dp 2 + dz 2 )+p 2 d^ 
A, PP + P~ X \ P + \ >zz = 



P fa - 



(1) 

(2) 
(3) 



Observe that (J2J is just the Laplace equation for A (in the Euclidean space). 
The 7-metric is defined by |5| 



A = - In 



R 1 + R 2 - 2m 
R1 + R2 + 2m 



3 2m 



(fli + R 2 + 2m) (R t + R 2 - 2m) 



(4) 
(5) 



where 



R\ = p 2 + ( z - m f Rl = p z + ( z + m )\ ( 6 ) 

It is worth noticing that A, as given by (jlj), corresponds to the Newtonian 
potential of a line segment of mass density 7/2 and length 2m, symmetrically 
distributed along the z axis. The particular case 7 = 1, corresponds to the 
Schwarzschild metric. 

It will be useful to work in Erez- Rosen coordinates [IT], given by 

p 2 = (r 2 — 2mr) sin 2 6 z — (r — m) cos 9, (7) 

which yields the line element as jU] 

ds 2 = Fdt 2 - F~ l [Gdr 2 + Hd6 2 + (r 2 - 2mr) sin 2 6d^ 2 ) , (8) 
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where 



G = ( , f- 2mr , . , Y~\ (10) 
V r 2 — 2mr + wr sin 9 I 



and 



h = - — ^ (ID 



(r 2 — 2mr) 7 

(r 2 — 2mr + m 2 sin 2 

Now, it is easy to check that 7 = 1 corresponds to the Schwarzschild metric. 

The total mass of the source is HD] M = 7m, and the quadrupole 
moment is given by 

Q = 2 m 3 ^ _ 7 2) (12) 

So that 7 > 1 (7 < 1) corresponds to an oblate (prolate) spheroid. 
2.2 The interior space-time 

The metric within the matter distribution bounded by the surface 

r = r s (13) 

is given by 

9u = ; 27 

9rr = _ / 2(l-7) A 7 2 -2 E l-7 2 

gee = -r*f^^ 

9 W = -r 2 / 2(1 ~ 7) sin 2 (14) 

where /, A, S and $ are functions whose specific form depends on the model 
under consideration. 

The two cases to be considered here are reported in [7j, namely 

1. The modified constant density Schwarzschild solution 



$(r,0) = f + ^V(r)sm 2 e 

with 

V(r) = l + - (r s -r), 
rs 

and 

B 2 - 



8np ss ' 

where p ss denotes the energy density in the spherically symmetric 
(7 = 1) 

2. The modified Adler solution 

f{r) = A + Br 2 

Cr 2 



A(r) = 1 + 



(A + Wr 2 ) 2/3 



Cr 2 C 2 r 4 

£(r,0) = 1 + ^tt + iH sin2 # 

(A + 35r 2 ) 2/3 4 (A + 35r 2 ) 4/3 

S(r,0) = (A + Br 2 ) 2 + - W(r) ... sin 2 9 
V 7 4 (A + 35r 2 ) 4/3 



with 



and 



m 



-l 



A = 



, , . 1 (r s -r) 



5m 
2r E 



1/2 
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2m (l 




2/3 



c 



1/3 



(27) 



s 



Before closing this section, two remarks are in order: 

1. Since we are considering the source described in (JT3j) as an initial state, 
the time derivatives of functions /, A, £ and $ will be in principle 
different from zero. 

2. Junction (Darmois) conditions are satisfied at the boundary r = 
-see [7] for details. 

3 The energy momentum tensor 

In order to give physical meaning to the components of the energy momentum 
tensor in coordinates (t, r, 9, ip), we shall develop a procedure similar to 
that used by Bondi ^Hj in his study of non static spherically symmetric 
sources. Thus, we introduce purely local Minkowski coordinates (r, x, y, 
z) (alternatively one may introduce a tetrad field associated with locally 
Minkowskian observers) defined by 



dr 



Pdt 



(28) 



dx 



jl-7^-l+7 2 /2 S (l-7 2 )/2 rfr 



(29) 



dy 



r/ 7(7-D$(l-7 2 )/2^ 



(30) 



dz 



rsin(^)/ 1 ^dip. 



(31) 



Next, since we are assuming that our source does not dissipate energy, 
then the covariant components of the energy momentum tensor, as measured 



by a local Minkowskian and comoving with the fluid observer, will be 



T, 



( p \ 

Vxx Pxy 

Pyx Pyy 

V o o o Pzz ) 



(32) 



where p is the energy density and in general p xx ^ p yy ^ p zz and p xy = p yx . 
We may write (|32|) in the form 

T^u — [p + Pzz) f [/ ~PzzV^u "I" (Pzz — Pzz) + (Pyj/ — Pzz) l^li> J r^'Pxyk(plu) > 

(33) 

where 7/^ denotes the flat space-time metric and 

% = ( 1, 0, 0, ) (34) 



0, 1, 0, 



(35) 



t = ( o, 0, 1, 



(36) 



The components of the energy-momentum tensor T Ml/ in (t, r, 8, (p) coor- 
dinates are linked to (jHB^l by 



where the non vanishing components of are 

A[ = P 



(37) 



(38) 



A x = jl-7 A -l+ 7 2 /2 S (l-7 2 )/2 
^ = r /'7(7-l) < |,(l-7 2 )/2 



(39) 
(40) 



r sin(0)/ 
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and the Lorentz matrices are given by 



-Vwi L\=Ll=5) + 



(r - \)wiWj 



(41) 



(42) 



where 

w 2 =w 2 x + w 2 r=-— L=, (43) 

Vl - w 

and w x and denote, respectively, the velocity of a fluid element along the 
x and y (r and 0) directions, as measured by our local Minkowskian observer 
as defined by (J25Jl - (J3T|l . Observe that we are considering the case w z = 0, 
which means that the system preserves the reflection symmetry (besides the 
axial symmetry). 

Then, (pT7jl readily gives 

T tt = / 27 r 2 (p + VxxWi + PyyWl + 2p xy W x W^j (44) 

(Tw x p + p xx w x A x + PyyWyA + p xy [w x A + w y A x ]) (45) 
T te = -r/^-^r x 

iTWyp + P XX W X A + PyyWyAy + ^ [^A + W X Ay\) (46) 

(r 2 u> 2 p + p^A 2 + p ra A 2 + 2 Pxy AA x ) (47) 

T r9 = r /(7-l) 2 A -l+7 2 /2,| ) (l-7 2 )/2 S (l-7 2 )/2 x 

{v 2 w x w y p + A [p xx A x + PyyAy] + p xy A 2 + A x A y ) (48) 

(r 2 ^^ + PxX A 2 + P yyA 2 y + 2p X yAAy) (49) 



with 



T w = r 2 f 2 ^ sin 2 (6)p zz , (50) 



A = (r ' 1} 2 W ^ , (51) 

A. 3 1 + (52) 
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A, - l + £z*S. (53, 



So that 



Tfiu — (p + Pzz) UfjJJp Pzz9fiv~\~ (Pxx Pzz) k^k u -\- (Pyy Pzz) I ~\~2pxyk(jJ>v) ) 

(54) 

where U^, and / M are obtained after applying the boost velocity and 

flEhgU to (jSl-iSl), 



Up = r ( /?, - Wx / 1 ^A- 1 +^ 2 / 2 S( 1 -^ 2 )/ 2 , - Wy rF<r'- 1 W 1 -r , V 2 , ) 

fc M = ( / 1 -7A- 1 +^ 2 /2 E (l-7 2 )/2A !r) r/ 7(7-l)$(l-7 2 )/2 A; Q ) 

^ = ( -Tw y f\ /l-7 A -l+7 2 /2 S (l-7 2 )/2 A; r /7(7-l)$(l-7 2 )/2 A ^ Q ) 



(55) 
(56) 
(57) 

In the static case w x = w y = A = and r = A x = A y = 1. The same 
holds after the system departs from equilibrium, on a time scale of the order 
of (or smaller than) the hidrostatic time scale (see Section 2.4). 



4 The Tolman mass and departure from equi- 
librium 

We shall now find an expression for the Tolman mass of the source described 
in the previous section, an then we shall evaluate it just after its departure 
from equilibrium, on a time scale of the order of (or smaller than) the hy- 
drostatic time scale. Also, since we are interested in the effects produced by 
small deviations from spherical symmetry, it will be assumed further that 
the source is quasi-spherical, and therefore we shall put 7 = 1 + e and will 
neglect terms of order 0(e 2 ) or smaller. 



4.1 The Tolman mass 

The Tolman mass is given by [Hj 
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+ 



JO 



2 77 



- 

dt 



-99 



a/3 



Tl - T 2 - Tl)drd6d(j) 



drd8d(j) (58) 



where C denotes the usual gravitational lagrangian density (eq.(10) in [Ej). 
Although Tolman's formula was introduced as a measure of the total energy 
of the system, with no commitment to its localization, we shall define the 
mass for any value of r, smaller or equal to rs, as 



rriT 



JO JO 

1 

*7T 



-g(T° - Tl - Tl - Tl)drdQd$ 



-99 



ai3 



d_ 

Of 



dC 



d(d(y/=gg"P)/dt) 



drd9d(j). (59) 



This extension of the global concept of energy to a local level [THj is suggested 
by the conspicuous role played by rar as the "active gravitational mass" . 

Indeed, it can be easily shown that the gravitational acceleration a 
of a test particle, instantaneously at rest in a static gravitational field, as 
measured with standard rods and coordinate clock is given by 



ttit 



(60) 



A similar conclusion can be obtained by inspection of the equation of hydro- 
static equlibrium (TOV) (in the static or quasi-static, spherically symmetric 
case case) [TT] . 

Even though these properties of Tolman's definition of mass are only valid 
in the spherically symmetric and static case, it is reasonable to assume that 
for small deviations from these conditions, the same role of active gravita- 
tional mass may be assigned to our expression. 

Let us now evaluate expression (joTJj) . After some lengthy calculations one 
finds: 



rriT 



2tt / / r 2 sin 6f 
JO Jo 

2 „/i i , .2 , , .2 

y 

+p zz + 2p xy (T 2 u x u y + A(A X + Ay)) 



2 7 +2 A 7 2 /2-l S (l-7 2 )/2 < | ) (l-7 2 )/2 



I ->( I • ,:i • Ji ) v + A 2 X + A 2 ) + Pyy (T 2 uj 2 y + A 2 + A 2 ) 

drdO 
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+ I f' r r 2 sin f 7 2 -4 7+ 2 A7 72-l S (l-7 2 )/2$(l-7 2 )/2 

4 Jo Jo 



/ 



+2 7 (l-7 2 )(-7 2 + 2 7 -3) 



2(2 - 3 7 + 7 2 )t + (7 2 - 2)| + (1- 7 2 )J + (1- 7 2 )^ 



E 
E 



r 



+ 2(7 2 -2)(7 2 -4) 



'A 



1 



(l-7 2 )(7 2 + l) 



E 



+ ( 7 2 -2)(2 - 3 7 ) 



/A 



+ (1 - 7 2 )(2 - 3 7 )|| + 7(1 - 7 2 )(27 - 3)^ | 



l-7 2 )(7 2 - 2)|| 



drdQ 



(61) 



We shall now evaluate this last expression at the very moment the system 
leaves the equilibrium, and for small non-sphericity 



4.2 Departure from equilibrium 

Let us now consider that our source, once submitted to perturbations, departs 
from equilibrium without dissipation. We shall then evaluate the system after 
such departure, on a time scale such that w x and w y remain vanishingly small, 
whereas their time derivatives though small, will be different from zero. 
Thus, just after leaving the equilibrium, the following conditions hold 

w x = w y = w Xji = w Vt i ~ 0, (i = r, 8, if) (62) 

w x , w y 7^ (small) (63) 

where dots denote derivative with respect to t. 

From now on, unless otherwise stated, all equations are evaluated at the 
moment the system starts to deviate from equilibrium. 

Then from (|44 )) -([5U )) . we obtain using (J52j) 

T te = T tr = (64) 
which implies, because of field equations (see for details) 

A = / = E = <i> = (65) 
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where for simplicity we write for 0{u) (as we shall do hereafter). 

Obviously, spatial derivatives of the above quantities will be also vanish- 
ingly small on the time scale under consideration. 



Also, from evaluation of the t-component of the conservation law T t 



we get (see [3]). 

p = (66) 

Then, in the approximation above, the expression for the Tolman mass 
becomes 

m T = 27T f r r 2 sin ^ /{72 - 27+2 ) A ( 7 V2-l) S (l-7 2 )/ 2<|) (l-7 2 )/2 

Jo Jo 

x (T ° - T\ - Tf - Ti)drd9d(f> 
A Jo Jo 



)/2 



X 



2(7 2 - 3 7 + 2)^ + ( 7 2 -2)| + (l- 7 2 ) (| + ^ 



drdOdcf) 



(67) 



We shall now write 7 = 1 + e, then assuming that non-sphericity is small 
(quasi-spherical approximation), the expression above reads: 

m T = 2tt T r r 2 sin 0/A (£ - 1/2) £ _£ $ _£ x (T ° - - T 2 2 - T*)drddd<f) 
Jo Jo 

+ I r r r 2 sin ^ / -(i +2e)A (e-i/2) S - e(|) - £ 

4 Jo Jo 



A / / A S $ N 



Within our two approximations (quasi-spherical and (J62|) - (}65|) ) it follows 
from (Pjl - fHHjl that: 



T ° - Ti - T| - T 3 3 = p + + p TO + p z 



(69) 



Also within these approximations, the time derivative of the tr component 
of the Einstein equations G r t = — 87rT rt , yields 



1 A 

r A 



+ 2e 



$' $' A 1 /A 
2$ + 2$2A + r [ A 



5 _ Z + ±^ _ ± + 
E / 2$i 2$U $ 



(70) 



13 



Feeding back ()69|) and (J70|) into ()68|). one obtains after some rearrangements 

m T = 2vr f r r 2 sin tf/A^ 1 ^^^ x (p + P<RB + p ra + p zz )drdd 
Jo Jo 

+ 2tt f J\ 3 smer 2e A^E~ 2£ ^ e [u x (p + Pxx )+u yPxy ]drde 



3 

- t 
1 

- 4 e 



x 



r 2 sin 6f- {1+2e) A^ 1/2) ^- e ^- e -drd9 



^sin^T^A^ 1 / 2 ^-^- 



$' $' A 1 /A 



_±(f& 

2$^2$2A%Ia E / ^ 2$ i 2$ 1 / $ 



drdO 
(71) 



4.3 The spherically symmetric case 

Before entering into the discussion of (|71|h it is quite instructive to analyze 
the spherically symmmetric situation. In this case one easily obtains from 



m T = 4vr / r f x (p + p xx + p yy + p zz )dr 

J 

+ 4vr f r ^A- 1 [dj x (p + Pxx )} dr (72) 
Jo 

with pyy = p zz . Now, the first term in (f72|l correspond to the active gravita- 
tional mass of the core of radius r interior to rs before the system leaves the 
equilibrium. Let us now assume that the system leaves the equilibrium and 
starts to collapse (expand), then u x < (u x > 0) decreasing (increasing) the 
value of the Tolman mass with respect to its value at equilibrium. Thus the 
second term in (ff^j) tends to stabilize the system, opposing any departure 
from equilibrium (contraction or expansion). This result was already known 

EE] 

Let us now get back to the non-spherical case to see how small non- 
sphericities change the whole picture. 



4.4 The non— spherical case 

Let us now turn back to ([71)1 . to infer what happens in the non-spherical case 
(i.e. e ytz 0), even though e « 1. Then neglecting higher order terms on e, 
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and keeping only the leading terms, we get 



m T = 2tt j f W r 2 sin ^/A (e " 1/2) S" e $" £ x (p + p xx + p yy + p zz )drd6 
Jo Jo 

+ 2tt r r r 3 sin fl/^A^ 1 ^- 2 ^^ [u x {p + p xx )\ drd6 
Jo Jo 

(73) 

In order to extract more information from ()73|) it is necessary to specify 
the source under consideration. We shall use the two configurations men- 
tioned in section 2. In both cases, / vanishes at the central region for 

> 2m. Thus, in the Schwarzschild-like models we have /(0) = 0, if 
2m/ Vs = 8/9. Since we know that any spherically symmetric static config- 
uration with constant p ss and locally isotropic pressure should satisfy the 
constraint n = 2M/r-£ < 8/9 we may assume that the system leaves the 
equilibrium for values of n close to 8/9. Then, an inspection of (|73|) shows 
that (if e > 0), for values of n approaching 8/9, the second term incresases 
due to the factor f~ 2e . 

We may now use Einstein equations to elucidate that for negative values 
of e the system may be unphysical. If n takes values close to the limit 
allowed by the model (8/9 for Schwarzschild-type model and 4/5 for Adler- 
type model), the critical values (/ — > 0) appear close to r = 0. Outside of 
this region the system is basically composed by an spherical incompressible 
fluid plus a perturbation in e. Thus, the physical or unphysical character of 
the model is determined by its behaviour close to these two limits. 

The energy density, in the limit r — > 0, for the Schwarzschild-type model 
and Adler-type model is given by expresions see (j3] for details) 



m j2e 

27rr| 



3 / 1 

+ e ' 



2 \ 3a /i_2»2_i 



(74) 



and 



3m 



P 



47rr| (1 — 2m/r s ) 



5m/2rs \ 



2e-2/3 



(1 - m/r s ) 2/3 + 



1/3 \^/l-2m/r s ) r (l-5m/2r 2 ) 1/3 

(75) 

respectively. From these ones, it is easy to show that if e < 0, the energy den- 
sity becomes negative as the system approaches to n —>■ 8/9 (Schwarzschild 
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case) or n — > 4/5 (Adler case) and positive for e > 0. Therefore, in both 
cases, positive energy conditions impose e > 0, and for very compact objects 
(close to the limit allowed by the model) the "stabilizying"term in the expres- 
sion for the active gravitational mass substantially increases in the central 
region, close to, but before, the maximum allowed value of n. 



5 Conclusions 

We have seen so far that (as expected), for high gravitational fields, important 
differences appear in the way that the spherical and non-spherical systems 
depart from equilibrium. This conclusion being true even for small non- 
sphericity. 

In the two models considered above, a factor multiplying the stabilizying 
term in the expression for the active gravitational mass, brings out those dif- 
ferences. Although for the two examples considered here, we have considered 
e > 0, it is obvious that there exist models with e < 0, in which case the sta- 
bilizying term may become very small, leading to highly unstable situations. 
In fact, it is worth noticing that important differences between the two cases 
(e > and e < 0) appear also in the behaviour of the exterior 7-metric ([§], 

Thus, on the basis of presented results we may conclude that whatever 
the model and the sign of e would be, any source of Weyl metric would evolve 
quite differently from the corresponding spherical source, as critical values of 
n are considered. 

In particular, in the examples examined, the departure from equilibrium 
appears to be affected by a sharp modification in the stabilyzying term of 
the active gravitational mass, as n approaches its maximum allowed value. 
This change makes the system more stable, hindering its departure from 
equilibrium. 

Before concluding, the following remarks are in order: 

1. Observe that all results above are valid only on a time scale of the 
order of (or smaller than) hydrostatic time scale. Within this time 
scale, the general form of the metric and the components of the energy- 
momentum tensor, is the same as in the static case (of course in the 
former case time dependence of variables has to be taken into account 
when taking time derivatives of such components) 
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2. In [7] it is shown that the two interiors described in Section 2.2, match 
smoothly (on r^) with the 7 metric, in the static case. That this is also 
the case after the system departs from equilibrium (on a time scale of 
the order of, or smaller than, hydrostatic time scale) is apparent from 
and comments after equation (57) and in the point 1 above. Of 
course in the fully dynamic case (beyond the hydrostatic time scale), 
(J65j) does not longer hold, and matching conditions have to be estab- 
lished. But this is out of the scope of this paper. 
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